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CN| , Abstract 

O , We briefly review the physics of gate operations between quantum dot spin-qubits and 

analyze the dynamics of quantum entanglement in such processes. The indistinguishable 
r- , character of the electrons whose spins realize the qubits gives rise to further entanglement- 

s' • like quantum correlations that go beyond simple antisymmetrization effects. We also 

summarize further recent results concerning this type of quantum correlations of indistin- 
guishable particles. Finally we discuss decoherence properties of spin-qubits when coupled 
to surrounding nuclear spins in a semiconductor nanostructure 

-xt- : 

^ ! 1 Introduction 

£NJ ' Quantum entanglement is one of the most intriguing features of quantum mechanics. It is 

one of the decisive properties of the quantum world which distinguishes it from classical 
| physics. Moreover, in the last decade, entanglement has proved to be the key ingredient 

_h ■ of the emerging field of quantum information theory 1 . 

In the beginning of modern quantum theory, the notion of entanglement was first 
• noted by Einstein, Podolsky, and Rosen [11], and by Schrodinger [12]. While in those 

days quantum entanglement and its predicted physical consequences were (at least par- 
tially) considered as an unphysical property of the formalism (a "paradox"), the modern 
perspective on this issue is very different. Nowadays quantum entanglement is to be seen 
as an experimentally verified property of nature providing a resource for a vast variety of 
novel phenomena and concepts such as quantum computation, quantum cryptography, or 
quantum teleportation. 

While the basic notion of entanglement in pure quantum states of bipartite systems 
(Alice and Bob) is theoretically well understood, fundamental questions are open concern- 
ing entanglement in mixed states (described by a proper density matrix), or entanglement 
of more than two parties [8, 9, 10]. The most elementary example for entanglement in a 
pure quantum state is given by a spin singlet composed from two spin-i-objects (qubits) 
owned by A (lice) and B(ob), respectively, 

^ (I T>a ® I i> B - I ® I T>u) . (i) 
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*To appear in the Proceedings of the E. Fermi School, "Quantum Phenomena of Mesoscopic Systems", 9 - 
19 July 2002, Varenna, Italy. 

1 For general introductions to the field of quantum information theory see, e.g., [1, 2, 3, 4, 5, 6, 7], for recent 
reviews on the issue of quantum entanglement see [8, 9, 10]. 
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For such a state, the state of the combined system cannot be described by specifying the 
state of Alice's and Bob's qubit separately. It is a standard result of quantum information 
theory [1, 4] that this property does not depend on the basis chosen in Alice's or Bob's 
Hilbcrt space. As we shall see below, the entanglement of such a quantum state (quantified 
by an appropriate measure) is invariant under (independent) changes of basis in both 
spaces. 

Physically measurable consequences of quantum entanglement of the above kind arise 
typically (but not exclusively) in terms of two-body correlations between the subsystems. 
In this case the effects of entanglement can typically be cast in terms of so-called Bell in- 
equalities [13] whose violation manifests the presence of entanglement in a given quantum 
state. Using this formal approach the physical existence of quantum entanglement (as op- 
posed to classical correlations) has unambiguously been verified by Aspect and coworkers 
for the polarization state of photons[14]. Moreover, quantum entanglement is an essen- 
tial ingredient of algorithms for quantum computation [2, 3, 4], in particular for Shor's 
algorithm for decomposing large numbers into their prime factors [15]. This problem is 
intimately related to public key cryptography systems such as RSA encoding which is 
widely used in today's electronic communication. 

Among the many proposals for experimental realizations of quantum information pro- 
cessing solid state systems have the advantage of offering the perspective of scalability, i.e. 
the integration of a large number of quantum gates into a quantum computer once the sin- 
gle gates and qubits are established. In the recent years, several proposals for using spins of 
electrons and/or nuclei in solid state systems have been put forward, starting with the work 
of Ref. [16], see Refs. [17, 18, 19, 20, 21]. Specifically, in Ref. [16] it was proposed to use 
the spin of electrons residing in semiconductor quantum dots as qubits [22, 23, 24, 25, 26]. 
The central idea is to store the information in the spin-degree-of freedom of the electron, 
while the manipulation of the spin, like single spin rotations and spin-spin interactions, 
are achieved via electrical gates acting on the charge degree-of-freedom of the electron. 
In this contribution we briefly review essential aspects of this proposal and analyze the 
dynamics of quantum entanglement in gate operations between qubits of this type. During 
such processes, the indistinguishable character of the electrons leads to entanglement-like 
quantum correlations which require a description different from the usual entanglement 
between distinguishable parties (Alice, Bob, ...) in bipartite (or multipartite) systems. In 
such a case the proper statistics of the indistinguishable particles has to be taken into 
account [24]. 

This article is organized as follows. In section 2 we give a brief overview on elementary 
results on entanglement between distinguishable parties. In section 3 we briefly review 
essential aspects of the quantum dot spin qubit proposal of Ref. [16]. After some gen- 
eral remarks we summarize in subsection 3.2 recent results on decoherence properties of 
electron spins interacting with surrounding nuclear spins in semiconductor nanostructures 
[53, 54]. In subsection 3.3 we give a brief description of the low-energy physics of lat- 
eral quantum dot molecules consisting of two spin-qubits. This discussion will lead us 
to the question of "entanglement" (or entanglement-analogous quantum correlations) be- 
tween indistinguishable particles, which are in our case here electrons in the joint system 
of two (or more) quantum dots [24]. The research and literature on quantum entan- 
glement has so far fairly concentrated on distinguishable parties; the issue of indistin- 
guishable (or, in other words, identical) objects has been addressed only very recently 
[24, 27, 28, 29, 30, 31, 32, 33, 34, 35, 36, 37, 38]. We therefore give also a short overview 
on this new direction of work in quantum information theory. In section 5 we discuss results 
on gate operations between quantum dot spin qubits with the emphasis on the dynamics 
of quantum entanglement. Specifically we present results on swap processes (interchanging 
the state of two qubits) and the "square root of a swap" generating a maximally entangled 
state from a tensor product. We close with a summary in section 6. 
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2 Quantum Entanglement between distinguishable par- 
ties 



We now give an introduction to basic concepts of characterizing and quantifying entan- 
glement between distinguishable parties. We concentrate on pure states (i.e. elements of 
the joint Hilbert space) of bipartite systems. We then comment only briefly on the case of 
mixed states (described by a proper density operator), and entanglement in multipartite 
systems. 

One of the most prominent examples of an entangled state was already given in the 
previous section, namely a spin singlet built up from two qubits. More generally, if Alice 
and Bob own Hilbert spaces Ha and Hb with dimensions m and n, respectively, a state 
\ip) is called nonentangled if it can be written as a product state, 

|V>) = \a) A ® \p) B (2) 

with \oi)a € Ha, \P)b € Hb- Otherwise \tp) is entangled. The question arises whether a 
given state \ip), expressed in some arbitrary basis of the joint Hilbert space H = Ha®Hb, 
is entangled or not, i.e. whether there are states \o)a and \0)b fulfilling (2). Moreover, 
one would like to quantify the entanglement contained in a state vector. 

An important tool to investigate such questions for bipartite systems is the biorthogonal 
Schmidt decomposition [1]. It states that for any state vector \xjj) <E H there exist bases of 
Ha and Hb such that 

r 

foft> = 5> (l«i> ® (3) 

i=l 

with coefficients Zi ^ and the basis states fulfilling (a,|aj) = (bi\bj) = 5^. Thus, each 
vector in both bases for Ha and Hb enters at most only one product vector in the above 
expansion. As a usual convention, the phases of the basis vectors involved in (3) can be 
chosen such that all Z{ are positive. The expression (3) is an expansion of the state \%p) 
into a basis of product vectors \a) <g> \b) with a minimum number r of nonzero terms. This 
number ranges from 1 to min{m,n} and is called the Schmidt rank of \ip). 
With respect to arbitrary bases in Ha and Hb a given state vector reads 

|^)=^M o6 |a)®|6) (4) 

a, b 

with an m x n coefficient matrix M. Under unitary transformations U a and Ub in Ha 
and Hb, respectively, M transforms as 

M h-> M' = U A MUl (5) 

with Ub being the transpose of Ub- The fact that there are always bases in Ha and Hb 
providing a biorthogonal Schmidt decomposition of is equivalent to stating that there 
are matrices Ua and Ub such that the resulting matrix M' consists of a diagonal block 
with only nonnegative entries while the rest of the matrix contains only zeros. For the case 
of equal dimensions of Alice's and Bob's space, m = n, this is also a well-known theorem 
of matrix algebra [39] . 

Obviously, |V>) is nonentangled, i.e. a simple product state, if and only if its Schmidt 
rank is one. More generally, the Schmidt rank of a pure state can be viewed as a rough 
characterization for its entanglement. However, since the Schmidt rank is by construction 
a discrete quantity it does not provide a proper quantification of entanglement. Therefore 
finer entanglement measures are desirable. For the case of two distinguishable parties, a 
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useful measure of entanglement is the von Neumann-entropy of partial density matrices 
constructed from the pure-state density matrix p = \ip}(tp\ [40]: 

E (\i>}) = - tT A (pa log 2 pa) = -tr B [p B log 2 p B ) , (6) 

where the partial density matrices are obtained by tracing out one of the subsystems, 
Pa/b = ^ t b/aP- With the help of the biorthogonal Schmidt decomposition of \ip) one 
shows that both partial density matrices have the same spectrum and therefore the same 
entropy, as stated in Eq. (6). In particular, the Schmidt rank of ip) equals the algebraic 
rank of the partial density matrices. \ip) is nonentangled if and only if the partial density 
matrices of the pure state p = IV'XV'I are also P ure states, and \ip) is maximally entangled 
if its partial density matrices are "maximally mixed", i.e. if they have only one non-zero 
eigenvalue with a multiplicity of min{m,n}. 

It is important to observe that the entanglement measure (6) of a given state does 
not depend on the bases used in Alice's and Bob's Hilbcrt space to express this state. This 
is because the trace operations in the definition of E(\ip)) are invariant under a change of 
bases (performed, in general, independently in both spaces). Therefore, entanglement in 
bipartite systems is a basis-independent quantity. 

Thus, the problem of characterizing and quantifying quantum entanglement for pure 
states in bipartite systems can been seen as completely solved. Unfortunately, the situation 
is much less clear for mixed states [8, 9, 10], and for multipartite entanglement. The main 
obstacle in the latter issue is the fact that the biorthogonal Schmidt decomposition in 
bipartite systems does not have a true analogue in the multipartite case. 

3 Quantum entanglement with electron spins in quan- 
tum dots 

We will now illustrate the phenomenon of quantum entanglement on the example of a 
specific (possible) realization of a quantum information processing system [16]. The pro- 
posal to be discussed below deals with qubits realized by the spins of electrons residing 
on semiconductor quantum dots. As we shall see in this and the following section, the 
indistinguishable character of the electrons gives rise to quantum correlations which are 
beyond entanglement between distinguishable parties. 

An study analogous to the one to be described below [24] has been carried out re- 
cently for the case of indistinguishable bosons (instead of electrons, i.e. indistinguishable 
fermions) in Refs. [41]. There the quantum dynamics (and generation of entanglement) 
between bosonic atoms trapped in neighboring potential minima of an optical lattice was 
studied and numerically simulated. 

3.1 General remarks 

An array of coupled quantum dots, see figure 1, each dot containing a top most spin 1/2, 
was found to be a promising candidate for a scalable quantum computer [16] where the 
qubit is defined by the spin 1/2 on the dot. Quantum algorithms can then be implemented 
using local single-spin rotations and the exchange coupling between nearby spins, see figure 
1. This proposal is supported by experiments where, e.g., Coulomb blockade effects, [42] 
tunneling between neighboring dots, [43, 42] and magnetization [44] have been observed 
as well as the formation of a delocalized single-particle state in coupled dots [45] . 

The basic mechanism for two-qubit gate operations within quantum-dot spin-qubit 
proposal are gated (i.e. time-dependent) tunneling amplitudes between the dots. Then 
the two-qubit dynamics is driven by the effective exchange coupling between the electron 
spins on different dots arising from a finite tunneling amplitude. 
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In the following subsection we shall briefly address the issue of decohcrcncc of spin- 
qubits. In subsection 3.3 we will summarize properties of a specific physical descrip- 
tion of quantum dot spin-qubits which are necessary for our discussion here. For more 
comprehensive reviews of quantum computing with electron spins in quantum dots see 
Refs. [46, 47, 48]. 

3.2 Decoherence properties 

The electron spin is a natural candidate for a qubit since its spin state in a given direc- 
tion, | |) or | J.}, can be identified with the classical bits |0) and |1), while an arbitrary 
superposition a\ |) +/3| J.) defines a qubit. In principle, any quantum two-level system can 
be used to define a qubit. However, one must be able to control coherent superpositions 
of the basis states of the quantum computer, i.e. no transition from quantum to classical 
behavior should occur. Thus, the coupling of the environment to the qubit should be 
small, resulting in a sufficiently large decoherence time T 2 (the time over which the phase 
of a superposition of |0) and |1) is well-defined). Assuming weak spin-orbit effects, the 
spin decoherence time T 2 can be completely different from the charge decoherence time (a 
few nanoseconds), and in fact it is known [50] that T 2 can be orders of magnitude longer 
than nanoseconds. Time-resolved optical measurements were used to determine T 2 * , the 
decoherence time of an ensemble of spins, with T 2 * exceeding 100 ns in bulk GaAs [50]. 
More recently, the single spin relaxation time T\ (usually T\ > T 2 ) of a single quantum dot 
attached to leads was measured via transport to be longer than a few fjs [51], consistent 
with calculations [52]. 

A particularly relevant and, from a fundamental point of view, also interesting mech- 
anism for electron spin decoherence in quantum dots (or other semiconductor nanostrcu- 
tures) is the hyperfine interaction with surrounding nuclear spins, which are ubiquitous in 
many materials such as GaAs. Thus, the Hamiltonian reads 

H = §-Y i A i h (7) 

i 

where S is the electron spin and Ii are the surrounding nuclear spins (both taken here to be 
dimensionless) The coupling constants Ai have dimension of energy and are proportional 
to the square modulus of the electron wave function at the location of the i-th nucleus. 
Since the electron wave function is spatially varying across the dot (a typical typical form 
for such envelope wave functions is given by a Gaussian), the coupling is inhomogeneous, 
i.e. different nuclei couple with different strength to the electron spin. This system was 
investigated in detail very recently in Refs. [53, 54]. Generally, the quantum dynamics 
according to (7) lead to a decay of the electron spin as measured in terms of its expectation 
values (S) (assuming an initially fully polarized spin, | (*S'} | = 1/2). This spin decay 
is generally a source of decoherence, i.e. loss of quantum information encoded in the 
qubit. However, a main finding in Ref. [53] is an exact analytical solution for the quantum 
dynamics for the situation of an initially fully polarized nuclear spin system with the 
electron spin pointing opposite to it. This exact solution can be evaluated explicitly in 
the limit of a large number ./V of nuclear spins being involved. As a result, the decay 
of the electron spin is in this case only of order l/N and occurs after a time interval of 
order TiN/A, where A = Ai cx 10 _5 eV for GaAs. In a typical GaAs quantum dot N 
is of order 10 5 . Therefore, for a sufficiently large number of nuclear spins being involved, 
the magnitude of the electron spin decay is small, and large decoherence times can be 
obtained. 

The findings of Ref. [53] have been corroborated and extended by numerical work in 
Ref. [54] . These investigations revealed a striking dependence of the electron spin dynamics 
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on the character of the initial state of the nuclear spin system (having a general degree of 
polarization) : If the initial state of the total system is a simple tensor product of the states 
each spin (both electron and nuclei), the time evolution depends strongly on the initial 
nuclear spin state. However, this behavior is very different from the situation where the 
electron is initially still uncorrelated with the nuclei, but the nuclei are among themselves 
randomly correlated, i.e. in a random superposition of all tensor product basis states of 
the nuclear system. Then the time evolution is very reproducible, i.e. it does (almost) not 
depend on the specific randomly generated initial state of the nuclear system; examples of 
such observations are shown in figure 2. Moreover, the time evolution of such a randomly 
correlated initial state mimics closely the average over the dynamics of all possible tensor 
product initial states. This observation can be understood as a self-averaging property 
[54] . The above results are illustrated in figures 2 and 3 

Another observation from the numerical studies in Rcf. [54] is the that the spin decay 
occurs more slowly if the nuclear spins are initially in a tensor product (i.e. non-entangled) 
state. This result can be expected to be of quite general nature and independent of the 
specific model (7): Disentangling the environment suppresses decoherence. 

The decay of the central electron spin S{t) in (7) is generally accompanied with the 
generation of entanglement between the electron and the nuclei [54] . This entanglement 
can be measured conveniently by the von-Neumann entropy of the partial density matrix 
where either the electron or the environment has been traced out from the pure-state 
density matrix of the total state [40], cf. Eq. (6). Tracing out the nuclear system we have 




This matrix has eigenvalues A± = 1/2 ± and the measure of entanglement reads 

E = — A + logA + — A_logA_. Thus, the formation of expectation values |(S(t))| < 1/2 
is a manifestation of the entanglement between the electron spin and the nuclear spin 
system. The maximum entanglement, E = log 2, is achieved if the electron spin has 
decayed completely as measured by the expectation values of its components, (S(t)) = 0. 
The generation of quantum entanglement between the electron spin and the nuclear spin 
system signaled by a reduced value of (S(t)} is a main and crucial difference between 
the quantum system studied here and its classical 'counterpart' described by a system of 
Landau-Lifshitz equations. 

3.3 Modeling of the double quantum dot system 

Let us now consider a system of two laterally tunnel-coupled dots having one electron 
each. Using an appropriate model [22, 23, 24] theoretical calculations have demonstrated 
the possibility of performing two-qubit quantum gate operations in such a system by 
varying the tunnel barrier between the dots. An important point to observe here is the 
fact that the electrons whose spins realize the qubits are indistinguishable particles [24]. 
Differently from the usual scenario of distinguishable parties (Alice, Bob, ...) the proper 
quantum statistics has to be taken into account when a finite tunneling between the dots 
is present [24, 27, 28]. 

The Hamiltonian of the double dot system is given by H = T + C, where C denotes 
the Coulomb repulsion between the electrons, and T = ^2 i=1 2 hi is the single-particle part 
with 

fH = ^ (ft + -/(n)) 2 + V(n) . (9) 

The single-particle Hamiltonian hi describes electron dynamics confined to the xy-plane in 
a perpendicular magnetic field B. The effective mass m is a material-dependent parameter. 
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The coupling of the dots (which includes tunneling) is modeled by a quartic potential 



rauj 



V(x,,) = =? x '-*> )■ + ,> , (10) 



which separates into two harmonic wells of frequency ojq (one for each dot) in the limit 
2a 3> 2a , where a is half the distance between the dots and do = yjTi/moj a is the effective 
Bohr radius of a dot. 

Following Burkard et al. [22] we employ the Hund-Mulliken method of molecular 
orbits to describe the low-lying spectrum of our system. This approach concentrates on 
the lowest orbital states in each dot and is an extension of the Heitler-London method also 
discussed in [22]. The Hund-Mulliken approach accounts for the fact that both electrons 
can, in the presence of a finite tunneling amplitude, explore the entire system of the 
two dots. Therefore this approach is suited to investigate the issue of entanglement- 
analogous quantum correlations between indistinguishable particles. In particular, the 
Hund-Mulliken approach includes two-particle states with both electrons being on the 
same dot. This issue of double occupancies will be discussed in more detail below. 

In the usual symmetric gauge A = B{— y,x,0)/2 the Fock-Darwin ground state of a 
single dot with harmonic confinement centered around r = (±a, 0,0) reads 



/raw fmu ( 2 2\\ 



•exp ( ±\y^T ) » ( n ) 



where I b — \fn~cfeB is the magnetic length, and the frequency w is given by u> 2 = + uj 2 l 
where u>l = eB/2mc is the usual Larmor frequency. From these non-orthogonal single- 
particle states we construct the orthonormalized states \A) and \B) with wave functions 

^ = R L 2 (<P+a-9<P-a), (12) 

^/l-2Sg-g 2 



with S being the overlap between the states (11) and g = (1 — \f\ — S 2 ) jS. For appropriate 
values of system parameters such as the interdot distance and the external magnetic field, 
the overlap S becomes exponentially small [22]. In this limit an electron in one of the states 
\A), \B) is predominantly localized around r = (±a, 0, 0). In the following we consider 
this case and use these states as basis states to define qubits, i.e. qubits are realized by 
the spin state of an electron in either orbital \A), or orbital \B). 

An appropriate basis set for the six-dimensional two-particle Hilbcrt space is given 
(using standard notation) by the three spin singlets 

\Si) = ^(4 t 4i-4i4t)|o>, (i4) 

\s*) = ^(4 T 4i + 4 T 4i)|o>, (is) 

\s 3 ) - ^ (4 T 4x - 4 T 4i) |o> , (i6) 

and the triplet multiplct, 

IT' 1 ) = c+^+JO), (17) 
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\T°) = ^(4t4i + 4i4 t )|0), (is) 
IT 1 ) = c+ T c+ T |0). (19) 

The three triplet states are degenerate (typically we can ignore possible Zeeman splittings 
[22]) and have the common eigenvalue, 



where we have defined 
and 



e T = 2e + V- , (20) 
e = (A\h\A) = (B\h\B) (21) 



VT = (T a \C\T a ) , ^+ = <S 1 |C|5i). (22) 

An important further observation is that, as a consequence of inversion symmetry along 
the axis connecting the dots, the Hamiltonian does not have any non-zero matrix elements 
between the singlet state IS3) and other states. Hence, IS3) is, independently of the system 
parameters, an eigenstate. The eigenvalues of the triplet and IS3), however, do depend on 
system parameters. The Hamiltonian acting on the remaining space spanned by |Si) and 
l^) can be written as 

«-* + it* +v t -(•£/» < 23 » 

where 

t H = -{A\h\B)- l -{S 2 \C\S l ) (24) 

and 

U H = (S 2 \C\S 2 )-V+. (25) 

The nontrivial part of (23) is a simple Hubbard Hamiltonian and can be identified as the 
Hamiltonian of a pseudospin- i-object in a pseudomagnetic field having a component Ur in 
the ^-direction and in the i-direction of pseudospin space. [Note that this pseudospin 
is not related to the spin degree of freedom which provides the qubit!] The space spanned 
by I Si) and IS2) contains the ground state of the system. The basis states themselves are 
eigenstates only in the case of vanishing tunneling amplitude tn where |Si) is the ground 
state. In all other cases, the ground state has an admixture of doubly occupied states 
contained in IS2). The energy gap between the triplet and the singlet ground state is 

st - sso = V + - V- - ^ + l - yJU* + ml . (26) 

A key challenge is the construction of systems composed of two coupled quantum dots 

which can be coupled to perform swap operations Usw, i-e. unitary two-qubit operations 

which interchange the spin states (qubits) of the electrons on the two dots By combining 
1/2 

the "square root" U sw of such a swap with other isolated-qubit manipulations one can 
construct a quantum XOR gate. A quantum XOR gate, along with isolated-qubit opera- 
tions, has been shown to be sufficient for the implementation of any quantum algorithm 
[55] . Hence a practical and reliable realization of a swap gate would be an important step 
toward the fabrication of a solid state quantum computer. A swap operation in the present 
system is a unitary transformation which turns a state having the qubits in different states, 
say, 

C + T 4j0) = i=(|T o ) + |Si)) , (27) 
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into a state where the contents of the qubits is interchanged, 

c+ l c+,\0) = ^=(\T°)-\S 1 )) . (28) 

These two states are eigenstates in the case V+ = V- and tn = for which the singlet- 
triplet splitting vanishes. 

As discussed in references [16, 22, 24], swapping may be achieved by the action of a gate 
that lowers the potential barrier between the quantum dots. This leads to exponentially 
larger values for both V + — V- and tn- It is adequate for our purposes to consider a model 
where V + = V- (consistent with the above limit of small overlap S), and the singlet-triplet 
splitting results entirely from t H . If the duration and amplitude of a tunneling pulse is 
adjusted appropriately, the relative phase between the singlet and the triplet state involved 
picks up a shift of 7r, and a swapping operation is performed. 

A finite tunneling amplitude necessarily leads to a finite probability for double occupan- 
cies of qubit states where both electrons are on the same dot. If double occupancy errors 
occur to any sizable extent as a result of the swapping process, any quantum computation 
based on this hardware is likely to fail. However, if the double occupancies are sufficiently 
rare after the swapping process, errors in the quantum computation can likely be corrected 
dynamically. An important observation is that the double-occupancy probability after the 
swap vanishes in the adiabatic limit, i.e. if the ramp time r of the quantum gate is such 
that Tl/t is much larger than the pseudospin splitting \JUfj + 16t 2 H . This follows since the 
non- adiabatic effects can arise only from the states \S\) and (S^), which have a non-trivial 
time evolution when the tunneling amplitude tn is time-dependent. Thus, the question 
of whether double occupancies are problematic for swap operations in the present system 
is reduced to the question of how close the motion of a spin-^-object in a time-dependent 
magnetic field is to its adiabatic limit. This will be investigated further in Sec. 5. 

The reduction of the dynamics to the time evolution of a two-level-system relies on 
the fact that the system has inversion symmetry along the x-axis in real space connecting 
the dots. This symmetry can be broken if odd powers of the particle coordinates Xi are 
added to the Hamiltonian (9) like for example the potential of a homogeneous electric 
field. However, the only additional matrix clement due to such terms in the Hamiltonian 
occurs in the subspace of double occupied states between the singlets IS2) and IS3). Thus, 
in the presence of an electric field £ = —eEj^i x i the Hamitonian acting on the singlet 
subspace spanned by \S\), IS2), \S$) reads 

H - 2e+^U H + V+ 

2t H 

-U H /2 F I (29) 





F 

with the real matrix element F = (5*2 |f IS3) and 

2X = (S 2 \C\S 2 ) - (S 3 \C\S 3 ) = 2(A\(A\C\B)\B) . (30) 

With a finite matrix element F the dynamics of the system is slightly more complicated, but 
also in this case the only coupling of the two-qubit states (27) and (28) to the subspace of 
double occupied states is provided by the tunneling amplitude tn- Therefore, with respect 
to the adiabaticity of the swapping process, the situation can be expected to be not very 
different from the one with inversion symmetry between the dots. This will be verified in 
Sec. 5. 

So far we have not considered a possible Zeeman coupling to the electron spin. This 
would not change the situation essentially since all states involved in the swapping process 
(|T°), |Si), IS2), and possibly IS3)) have the total spin quantum number S z = 0. 
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In the following section we give an elementary introduction to the theory of "entanglement- 
like" quantum correlations in systems of indistinguishable particles. We concentrate on 
the fermionic case and illustrate our findings on the above example of coupled quantum 
dots. 

4 Quantum Correlations between indistinguishable par- 
ticles 

For indistinguishable particles a pure quantum state must be formulated in terms of Slater 
determinants or Slater permanents for fermions and bosons, respectively. Generically, a 
Slater determinant contains correlations due to the exchange statistics of the indistinguish- 
able fermions. As the simplest possible example consider a wave function of two (spinless) 
fermions, 

*(n,f 2 ) = -j= [<Kri)x(r 2 ) - ^{r 2 )x(ri)] (31) 

with two orthonormalized single-particle wave functions <j){r), x(r). Operator matrix ele- 
ments between such single Slater determinants contain terms due to the antisymmctriza- 
tion of coordinates ("exchange contributions" in the language of Hartree-Fock theory). 
However, if the moduli of 4>(r), x{^) have only vanishingly small overlap, these exchange 
correlations will also tend to zero for any physically meaningful operator. This situation 
is generically realized if the supports of the single-particle wave functions are essentially 
centered around locations being sufficiently apart from each other, or the particles are sep- 
arated by a sufficiently large energy barrier. In this case the antisymmetrization present 
in Eq. (31) has no physical effect. 

Such observations clearly justify the treatment of indistinguishable particles separated 
by macroscopic distances as effectively distinguishable objects. So far, research in Quan- 
tum Information Theory has concentrated on this case, where the exchange statistics of 
particles forming quantum registers could be neglected, or was not specified at all. 

The situation is different if the particles constituting, say, qubits are close together 
and possibly coupled in some computational process. This the case for all proposals of 
quantum information processing based on quantum dots technology [16, 22, 23, 24, 25, 26]. 
Here qubits are realized by the spins of electrons living in a system of quantum dots. The 
electrons have the possibility of tunneling eventually from one dot to the other with a 
probability which can be modified by varying external parameters such as gate voltages 
and magnetic field. In such a situation the fermionic statistics of electrons is clearly 
essential. 

Additional correlations in many-fcrmion-systcms arise if more than one Slater deter- 
minant is involved, i.e. if there is no single-particle basis such that a given state of N 
indistinguishable fermions can be represented as an elementary Slater determinant (i.e. 
fully antisymmetric combination of N orthogonal single-particle states). These correla- 
tions are the analog of quantum entanglement in separated systems and are essential for 
quantum information processing in non-separated systems. 

As an example consider a "swap" process, as discussed in the previous subsection, exchanging 
the spin states of electrons on coupled quantum dots by gating the tunneling amplitude 
between them [22, 24]. Before the gate is turned on, the two electrons in the neighbor- 
ing quantum dots are in a state represented by a simple Slater determinant, and can be 
regarded as distinguishable since they are separated by a large energy barrier. When the 
barrier is lowered, more complex correlations between the electrons due to the dynamics 
arise. Interestingly, as shown in Refs. [22, 24], during such a process the system must 
necessarily enter a highly correlated state that cannot be represented by a single Slater 
determinant. The final state of the gate operation, however, is, similarly as the initial one, 
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essentially given by a single Slater determinant. Moreover, by adjusting the gating time 
appropriately one can also perform a "square root of a swap" which turns a single Slater 
determinant into a "maximally" correlated state in much the same way [24]. Illustrative 
details of these processes will be given in section 5. In the end of such a process the elec- 
trons can again be viewed as effectively distinguishable, but are in a maximally entangled 
state in the usual sense of distinguishable separated particles. In this sense the highly 
correlated intermediate state can be viewed as a resource for the production of entangled 
states. 

In the following we give an elementary introduction to recent results in the theory of 
quantum correlations in systems of indistinguishable particles [24, 27, 28]. This issue has 
attracted recently considerable interest and intense discussions [29, 30, 31, 32, 33, 34, 35, 
36, 37, 38]. The correlations to be discussed below are analogues of entanglement between 
distinguishable parties. However, to avoid confusion with the existing literature and in ac- 
cordance with Rcfs. [27, 28, 29], we shall reserve in the following the term "entanglement" 
for separated systems and characterize the analogous quantum correlation phenomenon in 
nonseparated systems in terms of the Slater rank and the correlation measure to be defined 
below. The results and concepts presented below and in the references given are in various 
respect complementary to a recent approach by Zanardi [31, 32] stressing the dependence 
of the notion of entanglement of the way a certain system is split into subsystems. 

For the purposes of this article we shall concentrate on elementary results for the case 
of pure states of two identical fermions. Results for mixed states and more than two 
fermions can be found in [24, 28] . Results for the case of identical bosons can be found in 
[30, 29, 28] 

We consider the case of two identical fermions sharing an n-dimensional single-particle 
space Jin resulting in a total Hilbert space A(H n <8> H„) with A denoting the antisym- 
metrization operator. A general state vector can be written as 

n 

H = E (32) 

a, 6=1 

with fermionic creation operators /+ acting on the vacuum |0). The antisymmetric coef- 
ficient matrix w a b fulfills the normalization condition 

tr(iBw) = ~, (33) 

where the bar stands for complex conjugation. Under a unitary transformation of the 
single-particle space, 

f+~Uf+U+ = U ha f+ , (34) 

w transforms as 

w i — ► UwU T , (35) 

where U T is the transpose (not the adjoint) of U . For any complex antisymmetric matrix 
n x n matrix w there is a unitary transformation U such that w' = UwU T has nonzero 
entries only in 2 x 2 blocks along the diagonal [27, 39]. That is, 



diag [Zi, . . . ,Z r , Zq] with Zi 



z, 
-Zi 



(36) 



Zi ^ for i e {1, . . . , r}, and Zq being the (n— 2r) x {n— 2r) null matrix. Each 2x2 block Zi 
corresponds to an elementary Slater determinant in the state \w'). Such elementary Slater 
determinants are the analogues of product states in systems consisting of distinguishable 
parties. Thus, when expressed in such a basis, the state \w) is s sum of elementary 
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Slater determinants where each single-particle basis state enters not more than one term. 
This property is analogous to the biorthogonality of the Schmidt decomposition discussed 
above. The matrix (36) represents an expansion of |u>) into a basis of elementary Slater 
determinants with a minimum number r of non- vanishing terms. This number is analogous 
to the Schmidt rank for the distinguishable case. Therefore we shall call it the (fermionic) 
Slater rank of \w) [27], and an expansion of the above form a Slater decomposition oi\w). 

We now turn to the case of two fermions in a four-dimensional single-particle space. 
This case is realized in a system of two coupled quantum dots hosting in total two electrons 
which are restricted to the lowest orbital state on each dot. In such a system, a simple 
correlation measure can be defined as follows [24, 27, 28]: For a given state (32) with a 
coefficient matrix to ab one defines a dual state \oJ a b) characterized by the dual matrix 

i 4 

c,d=l 

with e ahcd being the usual totally antisymmetric unit tensor. Then the correlation measure 
T)(\w)) can be defined as 2 



r)(\w}) = \(w\w)\ 



\ 

,b 



^ £ W ab W cd 



= |8 (w 12 W 34: + W13W42 + U>14?«23)| • (38) 



Obviously, rj(\w)) ranges from zero to one. Importantly it vanishes if and only if the state 
\w) has the fermionic Slater rank one, i.e. rj(\w}) is an elementary Slater determinant. 
This statement was proved first in Ref. [24]; an alternative proof can be given using the 
Slater decomposition of \w) and observing that 

detw = (l(w\w)] . (39) 



This is just a special case of a general relation expressing the determinant of an antisym- 
metric (2K) x (2K) matrix w by the square of its Pfaffian, 

/ 1 2K \2 

detw= \2K~Kl ^ £ ll '" l2K m ll2 ---w l2K _ ll2K J . (40) 

ii,...,i2Rr=l 

The quantity r](\w)) measures quantum correlation contained in the two-fermion state 
\w) which are beyond simple antisymmctrization effects. This correlation measure in under 
many aspects analogous to the entanglement measure "concurrence" used in systems of two 
distinguishable qubits [56]. These analogies are discussed in detail in [28] including also 
the case of indistinguishable bosons. An important difference between just two qubits, i.e. 
two distinguishable two-level systems, and the present case of two electrons in a two-dot 
system is that in latter system both electrons can possibly occupy the same dot while the 
other is empty. Therefore the total Hilbcrt space is larger than in the two-qubit system, 
and a generalized correlation measure becomes necessary. Furthermore, similar as in the 
two-qubit case, the correlation measure 77 defined here for pure states of two fermions has 
a natural extension to mixed fermionic and bosonic states [27, 28]. 

A convenient choice to make contact between the general state labels a, b, ... £ {1, 2, 3, 4} 
used here and the basis states of the preceding section is given by (1,2,3,4) = (A^,A I 



2 The notation n for the correlation measure was first used in Refs. [24, 27]. In Ref. [28] the notation C was 
used in order to stress the analogy to the "concurrence " introduced by Wootters [56] . 
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,B1,B I). With this convention, a state vector spanned by \S 2 ) and IS3) only has Wi 2 
and W34 as its only independent non-zero coefficients in w. Such a state lies fully in the 
subspace of double occupancies, and its entanglement is purely due to the orbital degrees 
of freedom: 

r) orb = 8|wi 2 W 3 4| • (41) 

On the other hand, a state spanned by |Si) and \T°) has no double occupancies and is 
entangled purely with respect to the spin degrees of freedom, 

Vspin = 8|wi 4 W23| • (42) 

For a general state vector, both kinds of correlations (orbital and spin) contribute to r/(w). 

A direct connection of the degree of quantum correlation r/ in the singlet ground state 
of the double dot system and its tunneling amplitude was found recently in [46]. The 
ground state reads 

K) = . * = ((1 + 0)1^) + (1 - 0)|S 2 » (43) 
^2(1 + (j) 2 ) 

with an energy given in (26). Here the "interaction parameter" 



has been introduced, and the correlation measure of the ground state is given by 

At zero tunneling tn = the ground state singlet does not contain any double occupancies 
and is degenerate with the triplet. At large tunneling tu 3> Uu both electrons are in 
symmetric orbital states differing in spin. This state is a single Slater determinant, and 
the correlation measure is zero. 



5 Dynamics of entanglement in quantum gate opera- 
tions 

We now continue with our investigation of the dynamics of the double quantum dot qubit 
swapping process generated by a time-dependent tunneling amplitude. 

Let us first consider the case of inversion symmetry along the axis connecting the dots. 
As explained in Sec. 3.3 this problem can be reduced essentially to the time evolution 
of a pscudospin-i -object in a magnetic field having a time-dependent component in the 
a;-direction of the pseudospin space. In the course of swapping, the triplet contribution to 
the incoming state (27) will just pick up a phase factor according to its constant eigenvalue, 
while the singlet contribution will mix with the other singlet IS2). Therefore, a finite prob- 
ability for double occupancies will necessarily arise during the swapping process. However, 
if these amplitudes can be suppressed sufficiently when the swapping is complete (as in 
the adiabatic limit), errors in the quantum computation can be avoided. Thus we are left 
with the question of how close the dynamics of our formal spin-l-object is to its adiabatic 
limit. We note that, in the adiabatic limit, no Berry phase occurs in the time evolution of 
the singlet states, since the motion of the formal spin is restricted to a plane. Hence the 
solid angle encircled in a round trip is strictly zero. 
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The integration of the Schrodingcr equation for our time-dependent two-level problem 
is in general non-elementary. However, there is a considerable body of literature, starting 
with early work by Landau [57], Zener [58], and Rosen and Zener [59], where particu- 
lar cases of this problem were reduced to well-known differential equation of mathematical 
physics such as the hypergeometric equation. This work was reviewed and generalized very 
recently in [60]. However, such an approach still works only for special time-dependent 
Hamiltonians, i.e., in the present context, only for special shapes of the tunneling pulse 
tii{t), and many quantities of interest are given by complicated non-elementary expressions 
which require numerical evaluation. For this reasons, and for the sake of brevity of our 
presentation, we shall resort to numerical integrations of the Schrodingcr equation. From 
such studies we will see that the range of adiabaticity is remarkably large. Our numeri- 
cal findings will be corroborated and made physically plausible by well-known analytical 
results for Landau-Zener-type transitions in simplified cases. 

To be specific, we consider a time-dependent tunneling of the form 

*"(*)= ., , cfU/ — ■ ( 46 ) 

1 + cosh(T/(2T)) 

This is a tunneling pulse which is switched on and off exponentially with a characteristic 
time r. It has a duration of T and an amplitude given by A (for T ^> r). Therefore 
this form is flexible enough to describe the essential features of a pulse. The exponential 
switching is motivated by the exponential-like dependence of the tunneling matrix element 
on external parameters [22] . 

A typical situation is shown in figure 4 for a switching time of r = AH/Uh, an amplitude 
of A = Uh/8 and the duration T adjusted to enable single swap operation. The figure 
shows the results of a numerical integration of the time-dependent Schrodingcr equation 
using the fourth order Runge-Kutta scheme. The time-dependent tunneling amplitude 
tii{t) is plotted (in units Uh) as a dotted line. The square amplitude of the incoming 
state (27) and the outgoing state (28) are shown as thick lines. The square amplitudes 
of the singlets |<Si) and IS2) are denoted by \<pi\ 2 and \(f2\ 2 , respectively, and plotted as 
long-dashed lines. The probability of double occupancies is given by | 1 2 - As one can see 
from the figure, this quanity is finite during the swapping process but strongly suppressed 
afterwards. The measure of entanglement rj(t) is also shown in the figure. It is zero for 
the non-entangled incoming and outgoing state, and achieves its maximum value of almost 
unity in the middle of the process. This quantifies and shows explicitly the entanglement 
of the quantum state during the swapping process. 

The probability \ip2\ 2 for double occupancy after switching off the tunneling depends 
on the switching time r, the amplitude A and also on the duration T of the tunneling 
pulse, i.e. on the exact time when the switching off sets in. However, our numerics suggest 
that there is an upper bound for \(fi2\ 2 at given r and A. In the above example the double 
occupancy probability after the swapping process is smaller than 10 -10 , which is a very 
tiny value. A typical order of magnitude for the double occupancy probability is 1CT 6 for 
amplitudes A < Uh and switching times r > 4Ti/Uh- In fact, also larger values of A (being 
still comparable with Uh) can be possible, leading to double occupancy probabilities of 
the same order, while this probability significantly increases if r becomes smaller than 
4Ti/Uh- Thus, this value characterizes the region where the motion of the system is close 
to its adiabatic limit and is remarkably small on the natural time scale of the system 
given by H/Uh, while adiabatic behavior is in general expected for a particularly slow 
time evolution. 

This large range of quasi-adiabatic behavior can be understood qualitatively by con- 
sidering a simplified situation where the tunneling is switched on and off linearly in time 
and is constant otherwise. Then, non-adiabatic effects can occur only during the sharply 
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defined switching processes. For simplicity, we consider the first switching process only 
where the tunneling has the time dependence t H = (A/r)i, t e [0, r]. To enable analytical 
progress let us further assume t G [—00,00], which should lead to an upper bound for the 
probability of non-adiabatic transitions due to the switching. This problem was consid- 
ered a long time ago by Landau [57] and by Zener [58] . The result of reference [58] for the 
probability of non-adiabatic transitions reads 

P„ad = e~ a (47) 

with an adiabaticity parameter 

We see that the probability for non-adiabatic transitions is exponentially suppressed with 
increasing switching time r. This exponential dependence explains qualitatively the above 
observation of a large range of quasi-adiabatic behavior. To obtain an estimate for a 
nonlinear switching one may replace the ratio (A/r) in the denominator of (48) by the 
maximum time derivative of the tunneling t H (t) (giving a = irUfj / '37i(A / V) for the pulse 
(46)). 

A similar exponential dependence of the probability for non-adiabatic transitions on 
the switching time r was also found analytically by Rosen and Zener [59] for a particular 
two-parametric pulse of the form 

t H {t) = A/cosh(t/r). (49) 

In this case non-adiabatic transitions occur with a probability 

Pnad = sin 2 (Ar/(2ft)) / cosh 2 (U h t/(2H)) . (50) 

To illustrate the behavior in the strongly non-adiabatic case we have plotted in figure 5 
and |</?2| 2 for the same situation as in figure 4 but with a four times smaller ramp 
time of only r = H/Uh- In this case small oscillations occur in the time-evolution of 
these two quantities during the tunneling pulse, which can be understood in terms of the 
eigenspectrum at a given tunneling tn — A. Additionally, a sizable double occupancy 
probability of about 0.005 remains after the pulse, as shown in the inset. 

Figure 6 shows a square root of a swap, which is obtained from the situation of figure 
4 by halting the duration T of the tunneling pulse. The resulting state is a fully entangled 
complex linear combination of the states \S\) and IT ), or, equivalently, of the incoming 
state (27) and the outgoing state (28) of the full swap. Again, the weight of the doubly 
occupied state \S%) is strongly suppressed after the tunneling pulse. As a consequence, Eq. 
(41) implies that r) or b = after completion of switching, while 77 = t] sp i n — 8\wi4U>23\ = 1. 
This shows that the entanglement of the two electrons is entirely in the spin (and not in 
the orbital) degrees of freedom after switching. 

Let us finally consider swapping processes when the inversion symmetry along the axis 
connecting the dots is broken. Such processes are governed by the Hamiltonian (29) in the 
presence of a finite matrix element F. Our numerical results are in this case qualitatively 
the same as before with the admissible switching times r slightly growing with increasing 
F. In figure 7 we illustrate our findings for a comparatively large off-diagonal element 
F = OAUh- The additional Coulomb matrix element is X — 0.2Uh, and the parameters 
of the tunneling pulse are r = 8H/Uh and A = Uh/S with a duration T appropriate for 
a single swapping. As a result, a clean swapping operation can be performed also in the 
absence of inversion symmetry. 

We note that the Hund-Mulliken scheme used here is restricted to the low-energy sector 
where only the lowest single-particle energy levels (with typical spacings Se) are kept. For 
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this scheme to be valid also in a switching process, we need to require that time-dependent 
changes must be performed adiabatically also with respect to the time scale set by h/Se, 
i.e. we need r > h/Se[22]. On the other hand, to suppress double occupancy errors we 
have seen that the adiabaticity parameter a of Eq. (48) must exceed unity, implying that 
t > 8UA/(ttU h ). Thus, the adiabaticity condition for switching becomes more generally, 



There arc now two particular cases we can distinguish. First, if the effective Coulomb 
charging energy exceeds the level spacing, i.e. Uh > Se, we obtain r m i„ = h/Se, since for 
consistency we have A < Se. Thus, when the switching is adiabatic with respect to the 
scale set by Se, errors due to double occupancy are automatically excluded. In the second 
case with Uh < V ASe < Se ("ultrasmall quantum dots"), we obtain T min = &KA/(ttU h ), 
which means that the overall condition for adiabaticity is determined by the no-double 
occupancy criterion. Using typical material parameters for GaAs quantum dots [43], we 
can estimate[22] that T min is of the order of 50 ps. 

6 Summary 

We have given an overlook on the physics of gate operations between quantum dot spin- 
qubits emphasizing the dynamics of quantum entanglement between such objects. For 
these purposes, the low-energy physics of a lateral pair of identical quantum dots can be 
effectively reduced to a (time-dependent) two-state problem, In particular we have ana- 
lyzed the swap and "square root of a swap" operations by means of numerical simulations 
of the time-dependent Schrodinger equation. One of the main conclusions is that double 
occupancy error will not be a fatal problem for the operation of such quantum gates. This 
is due to a surprisingly large range of the quasiadiabatic regime in the quantum dynam- 
ics of the aforementioned two-state system Our numerical results are corroborated and 
physically made plausible by analytical results on related Landau-Zener-type transitions. 
We have also reported on recent results on electron spin decoherence in semiconductor 
nanostructures due to hyperfinc interaction with surrounding nuclear spins. 

Finally, the indistinguishable character of the electrons involved in the realization of 
quantum dot spin-qubits leads to the interesting question of entanglement-analogous quan- 
tum correlations between indistinguishable particles. Here we have provided an introduc- 
tion to this issue, which has been addressed in the literature since only very recently and 
constitutes a new direction of work in quantum information theory. 
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high-g layer quantum well 

Figure 1: Quantum dot array, controlled by electrical gating. The electrodes (dark gray) 
define quantum dots (circles) by confining electrons. Additional gates and (dc and ac) external 
magnetic fields enable single-qubit rottaions. Two-qubit operations are achived by gating the 
effective tunneling amplitude between neighboring quantum dots. This introduces an effective 
time-dependent exchange coupling between the electron spins. 
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Figure 2: Upper panels: (S z (t)} for a system of size N = 19 being initially in a randomly 
correlated nuclear spin state in the subspace with J z = 7 where J = S+J2i h denotes the total 
spin of the system. The two panels represent two different randomly chosen initial conditions. 
Lower panels: Analogous data for N = 14 and a completely unpolarized nuclear spin system 
(J z = —1/2). In both cases the simulation data does practically not depend on the initial 
condition. 
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Figure 3: Upper panels: time evolution of the electron spin (S z (t)} for a system with 14 nuclear 
spins being initially in an uncorrelated tensor product state in the subspace J z = 9/2. The 
time evolution of the electron spin depends strongly on the initial state of the nuclear spins. 
Lower left panel: data of the same type as above but averaged over all possible uncorrelated 
initial states with J z = 9/2. Lower right panel: (S z (t)) for the same system being initially in 
a randomly chosen correlated state. This time evolution closely mimics the avergaged data of 
the lower left panel. 
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Figure 4: A swap process as a function of time. The tunneling amplitude tn(t) is plotted 
(in units of Uh) as a dotted line. The square amplitude of the incoming state (27) and 
the outgoing state (28) are shown as thick lines. The square amplitudes of the singlets \S\) 
and IS2) are denoted by |v?i| 2 and \(f2\ 2 , respectively, and plotted as long-dashed lines. The 
measure of entanglement rj(t) is also shown. 
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Figure 5: The square amplitudes of the singlet states and IS2) for the same situation as in 
figure 4, but with a four times smaller ramp time of only r = Ti/Uh- The inset shows \(f2(t)\ 2 
on a magnified scale. The dynamics of the system is clearly in the non-adiabatic regime. 




Figure 6: A square root of a swap, which is obtained from the situation of figure 4 by halfing 
the pulse duration T. The probability of double occupancies is again strongly suppressed after 
the tunneling pulse. The resulting state is a fully entangled complex linear combination of 
| Si) and \T°), or, equivalently, of the incoming state (27) and the outgoing state (28) of the 
full swap. The quantum mechanical weigths of the latter states are plotted as thick solid lines 
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Figure 7: A swapping processes in the absence of inversion symmetry along the axis connecting 
the dots. The square amplitudes of the singlet states \Si), i G {1,2,3}, are denoted by \<£i\ 2 ■ 
The additional matrix elements entering the Hamiltonian (29) are X = 0.2Uh and F = OAUh- 
The parameters of the tunneling pulse are r = 8H/Uh and A = Uh/8 with a duration T 
appropriate for a single swapping. As a result, a clean swapping operation can be performed 
also in the absence of inversion symmetry. 
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